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Abstract. We obtain a natural extension of the Vlasov-Poisson system for stellar 
dynamics to spaces of constant Gaussian curvature k ^ 0: the unit sphere S^, for 
K > 0, and the unit hyperbolic sphere for k < 0. These equations can be easily 
generalized to higher dimensions. When the particles move on a geodesic, the system 
reduces to a 1-dimensional problem that is more singular than the classical analogue 
of the Vlasov-Poisson system. In the analysis of this reduced model, we study the 
well-posedness of the problem and derive Penrose-type conditions for linear stability 
around homogeneous solutions in the sense of Landau damping. 


1. Introduction 

The Vlasov-Poisson system models the density change of galaxies in a cluster of 
galaxies, stars in a galaxy, or particles in plasma. The galaxies, stars, or particles are 
assumed to be identical to each other, while collisions, relativistic effects, and magnetic 
helds are neglected. If ignoring these effects is physically unreasonable in certain prob¬ 
lems, some related models can be used, such as the Vlasov-Maxwell [1] , Einstein-Vlasov 
[30], or Vlasov-Manev systems [5]. We focus here on the case of stellar dynamics, as¬ 
suming that collisions do not occur and relativistic effects can be neglected. Under these 
assumptions, the evolution of stars or galaxies is usually modelled in the framework of 
kinetic theory by the Vlasov-Poisson system in Euclidean space, given by the equations 
d d d 

^/(^,x,v) + V— /(t,x,v) + F(t,x) —/(t,x,v) = 0, 

F{t,x) = -[ — ^p{t,y)dy, p{t,x) := [ /(t,x,v)dv, 

where U represents the potential. For an initial value problem, a value /(0,x, v) is 
given to determine the distribution function f(t, x, v) of the celestial objects at position 
X G with velocity v G at time t G M. Due to the dependence of the potential on 
the density p, the system reduces to a non-linear partial differential equation, which is 
difficult to understand. Even the global existence of solutions remained an open problem 
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for decades. The first to attack it was R. Kurth, who showed the local existence of 
solutions in 1952 [2T]. In 1977, J. Batt proved global existence for spherically symmetric 
solutions [1] and, in 1985, C. Bardos and P. Degond completed the global existence of 
solutions under the assumption of small initial data [3]. Finally, the existence of global 
solutions with general initial data was achieved by K. Pfaffelmoser in 1990 [28] and, 
independently, by P.L. Lions and B. Perthame in 1991 [23]. Recently, C. Mouhot and 
C. Villani showed that these equations exhibit Landau damping, a remarkable result 
that inspired a direction of research we took here, [26] . 

In this paper we broaden the scope of the Vlasov-Poisson system to spaces of non¬ 
zero constant Gaussian curvature within the framework of classical mechanics, without 
involving special or general relativity. On small scales, the curvature of the physical 
space is negligible, but we cannot exclude the possibility that the universe is hyperbolic 
or elliptic on the large scale. By extending the study of the Vlasov-Poisson system 
to spaces of non-zero constant curvature we could, on one hand, obtain a better un¬ 
derstanding of the flat case by viewing the Vlasov-Poisson system as the limit of its 
counterpart in curved space when the curvature tends to zero. On the other hand, we 
might be able to decide whether the universe is curved or not. Indeed, if a certain 
solution of the density function occurs only in, say, flat space but not in hyperbolic and 
elliptic space, and such behaviour is supported by astronomical evidence, then we could 
claim that space is Euclidean. Even if only for these two reasons alone, the extension 
of the Vlasov-Poisson system to spaces of constant curvature deserves a detailed study. 

The derivation of the classical system for stellar dynamics uses the Newtonian equa¬ 
tions of the gravitational iV-body problem. To generalize the Vlasov-Poisson system 
to elliptic and hyperbolic space, we employ a meaningful extension of the Newtonian 
V-body problem to spaces of constant Gaussian curvature. Although the idea of such 
an extension belonged to Bolyai and Lobachevsky in the 2-body case [21E], a suitable 
generalization of the classical Newtonian system was only recently obtained, mm 
Suitability is meant here in a mathematical sense, since there are no physical ways of 
testing the new system. More precisely, the potential of the “curved problem” in its 
most simple setting (the case of one body moving around a hxed attractive centre—the 
so-called Kepler problem), not only recovers the Newtonian case as the curvature tends 
to zero, but also inherits two properties of the classical potential: (i) it is a harmonic 
function, i.e. satishes Laplace’s equation; (ii) all of its bounded orbits are closed. More¬ 
over, the properties that have been so far discovered for the extension of the V-body 
problem to spaces of constant curvature ht nicely with what is known in the classical 
case POromsiillDE]. 

Our derivation of the Vlasov-Poisson system is based on these equations of the curved 
V-body problem and on Liouville’s theorem. We could also regard these equations in a 
mean held approximation to obtain, at least formally, the Vlasov equation in spaces of 
constant curvature. Both these approaches use coordinate systems (extrinsic or intrin¬ 
sic). An approach in the spirit of geometric mechanics to dehne the Vlasov equation, 
without the use of coordinates, is due to J. Marsden and A. Weinstein [25] . This method 
regards the Vlasov equation as a conservation of a distribution function along phase- 
space trajectories of a given Hamiltonian. When the potential in the Hamiltonian is 
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given by Poisson’s equation, we naturally obtain a Vlasov-Poisson system on symplec- 
tic manifolds and, in particular, on spaces of constant curvature. The analysis of such 
systems on manifolds appears to be difficult since the transport part does not seem to 
be well understood. For this reason, at the current stage of our research, we focus on 
the special case when the particles move along a geodesic of the manifold. As geodesics 
are invariants of the equations of motion, such conhgurations are also invariant. By im¬ 
posing this restriction, we derive a new (1-dimensional) reduced model given by system 
(4.2) in Proposition 4.1, which provides our hrst notable result. 

This model is different from the classical 1-dimensional Vlasov-Poisson system [2] 
and the Vlasov-Hamilton Mean Field equations [IB] because the interaction potential is 
more singular. However, the mathematical challenges are similar to those of the Vlasov 
equation given by the Dirac potential (also called gyrokinetic), derived in [29], [T8] . 
121 . which is important in the analysis of plasma fusion. Given the singularity in the 
potential, we focus (as in [19]) on the existence of strong solutions, which we construct 
in the class of Gevrey functions. This is our second notable result, and we state it in 
Theorem 15.11 

Following G. Villani [BT], we then derive Penrose-type conditions for our 1-dimensional 
Vlasov-Poisson system such that the homogeneous states are linearly stable, i.e. require 
convergence for density and force. The convergence (or decay) rate is exponential when 
the model is reduced from the sphere. We want to emphasize that this decay becomes 
only algebraic for the reduced hyperbolic system. This phenomenon occurs because of 
the low frequencies when we work on the whole real line. In this context, our hnal 
notable result, which states that the linearized Vlasov-Poisson system along a geodesic 
has stable solutions in the sense that they exhibit the Landau damping phenomenon, 

the former dealing with the sphere and the latter 


is stated in Theorems |6.1| and 6^ 
treating the hyperbolic case. 

Our paper is organized as follows. In Section 2, we lay the background for further 
developments. After introducing some notation, we provide a solution representation 
to Poisson’s equation based on the solution of the Laplace-Beltrami equation obtained 
in [7], for the sphere, and [8], for the hyperbolic sphere. In Section 3, we derive the 
Vlasov-Poisson system on 2-dimensional curved spaces (the unit 2-sphere and the unit 
hyperbolic 2-sphere), using the extension of the Newtonian equations mentioned above, 
and write this system in Hamiltonian form. The generalization of these equations to 
higher dimensions is straightforward. Then we obtain an expression of the gravitational 
held generated by the bodies (assumed to be point particles), which we apply to a 1- 
body example to illustrate the use of symmetry and the agreement between our original 
assumptions and our derived equations. For simplicity, in Section 4 we choose initial 
data such that the particles move on a geodesic of the sphere or hyperbolic sphere. 
Since the geodesics are invariant sets for the equations of motion, this approach allows 
us to use the symmetry of the conhguration to simplify the system. We continue under 
this assumption for the remainder of the paper, and in Section 5 we prove the existence 
of solutions to our Vlasov-Poisson system. Finally, in Section 6, we state and prove our 
stability results mentioned above. 
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2 . Background 


Since we are interested only in qualitative properties of solutions, we will choose from 
the start the physical units such that the gravitational constant is 1. In the discrete 
case, when dealing with a hnite number of point masses, the qualitative study of the 
motion can be reduced to the unit sphere, for positive curvature, and the unit hyperbolic 
sphere, for negative curvature, mm- Since this framework can be used without loss 
of generality in our qualitative studies as well, we will also employ it here. 

Let represent the unit 2-sphere, 

= {{xi,X2,X3) \xl + xl + xl = 1 }, 


for positive curvature, and the unit hyperbolic 2 -sphere (i.e. the upper sheet of the unit 
hyperboloid of two sheets), 

EI^ = {{xi,X 2 ,Xz) \ xl + xl- xl = -1, Xg > 0}, 


for negative curvature, where is embedded in the 3-dimensional Euclidean space 
and is embedded in the 3-dimensional Minkowski space in which all coordinates 
are spatial. The origin of the coordinate system lies at the centre of the sphere and the 
hyperbolic sphere. If k denotes the Gaussian curvature, the signum function is dehned 
as 


a = 



for K > 0 
for K < 0 . 


Following m and [ 12 ], we introduce unihed trigonometric functions. 


( 2 . 1 ) 


snx 


sinx, for k > 0 
sinhx, for k < 0 , 


csnx 


cosx, for K > 0 
coshx, for K < 0 , 


snx 

tnx := -, 

csnx 


csnx 

ctnx := -, 

snx 


in order to treat both and simultaneously. In this framework, we can write the 
geodesic distance between the points a = (oi, 02, Og) and b = ( 61 , b2, &g) of as 


( 2 . 2 ) 


d(a, b) = csn ^(cxa-b). 


where the operation ■ stands for the scalar product 


a ■ b — aibi 0,2^2 T 

which is inherited from the space in which the manifold is embedded, i.e. for but 
for We dehne the norm of the vector a as 

||a|| = \al + al + 

The gradient of a scalar function f = f(x), with x G M^, is given by 
(2.3) Vxf = (c^xif) Gi -|- {dx2^) 62 + {(jdxJ) eg. 
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where 61 , 62,63 denote the elements of the canonical base in Cartesian coordinates. 
The divergence of a vector function F = Fi(x) 6 i + F 2 (x )62 + F 3 (x) 63 , where x := 
TiOi + 0:262 + 0:363 G M^, has the form 

(2.4) divx F = d^^Fi + d^^F 2 + d^^F^. 

The Laplace-Beltrami operator on M^, applied to a scalar function f as dehned above, 
can then be written as 


(2.5) Axf — dx^xif + dx2X2f + ^dxgxsf- 

Notice that the operators Vx, diVx, and Ax act on functions dehned on the manifold 
after these functions are extended to the 3-dimensional ambient space by replacing the 
variable x with x/||x||. Once the computations are performed, the resulting functions 
can be restricted to the manifold again by imposing the condition ||x|| = 1 , which 
dehnes M^. 

Recall that the tangent space at x G is given by 

TxM^ = {v I X ■ V = 0 } 

and the tangent bundle TM^ is the union of all tangent spaces, i.e. 

TWf = IJ TxM^ = {(x,v) I X G x ■ v = O} . 

xeM2 

To parametrize M^, let /m 2 be the interval [ 0 , 7 r] if but [0, cx)) if 

Then the position of a particle at x G can be represented as 


( 2 . 6 ) X = a:i 6 i -|- 0:262 -1- 0:363 = (sn a cos 9)ei + (sn a sin 9 )e 2 + (csn q;) 63 , 

where 9 G [0, 27r) and a G /m 2 . We dehne 6^, 6^, and ee to be the orthogonal vectors 
that satisfy the equation 



0^ 


sn a cos 9 

sn a sin 9 

csn a 


ei 

(2.7) 

_e0_ 

— 

csn a cos 9 
— sin 9 

csn a sin 9 
cos 9 

—asncK 

0 


e2 

_e3_ 


When the parameters a and 9 depend on time, the velocity of a particle at x is 


V = OOa 6q, -(- {U 0 sn a) eg, 

where Ua ■= a and ug := 9, and the upper dot denotes the derivative relative to time. 
The acceleration is then given by 

a = —a (ca^ -|- sn ‘^a)er + {ua — ujjsna csn a)ea + (coesn a + 2 ua uje csn a)eg. 


The natural Riemannian or pseudo-Riemannian metric on has the form 


[dij] 




1 0 
0 sn ’ 


with inverse 



ij=l,2 


1 

0 


0 

1 


sn 2 q;-I 
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Then, in the local coordinates q := cte^ + 9eg, the gradient of f = f(a, 9) := f(x(q)) at 
q G is given by 

(2.8) 


V,f=(5J)e„+ -e», 


sn a 


the Laplace-Beltrami operator takes the form 


Aq f = diVq(Vqf) = ctn a + d^af H- deef, 

sri cx. 

where the divergence of a vector fnnction 

F = F(a, 9) := F(x(a, 9)) = Fa{a, 9)ea + Fg{a, 9)eo 
at q G is expressed as 

divq F = Fa ctn a + daFa H- dgFg. 

sna 


The volnme form on is 


= sn a dad9. 


3 . Equations of motion, gravitational potential and the curved 

Vlasov-Poisson system 

In this section we first derive the equations of motion for a particle moving on the 
manifold M^, then introduce the gravitational potential and the gravitational force 
function, and finally obtain the Vlasov-Poisson system on and Although, to fix 
the ideas, we work here only in the 2-dimensional case, our equations can be extended 
to any hnite dimension. 


3.1. Equations of motion. Recall that for (x, v) G TM^, / = /(f,x, v) denotes the 
phase space density, p = p(t, x) the spatial density, and U = U(t, x) the gravitational 
potential function. 


Proposition 3.1. In extrinsic coordinates having the origin at the centre o/M^, the 
eguations of motion for a particle with position x G and velocity v = uiOi - 1 - V2G2 + 
11303 G TM^, under the effect of a potential function V : M x —)■ M, are 


(3.1) 


X = V, 

V = V:x_U(t, x) — (t(v ■ v)x. 


Proof. The proof of this result can be found in Section 3.4 of m for a finite number 
of masses for the Newtonian potential function. We will adapt it here to the form of 
the potential function U given in (3.10) in the case of one body. We need to use hrst 
a variational result from constrained Lagrangian mechanics (see, e.g., HI). since the 
motion of the particle is restricted to M^. Under such circumstances, let L represent 
the Lagrangian and gfx) = 0 be the equation expressing the constraint, i.e. 


L(«,x,v) = r(x,v) - V^(«,x) = -cr(v ■ v)(x. x) + t;(«,x), 


(3.2) 
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and 

(3.3) g{x.) = xl + xl + axl — a = 0. 

The factor (j(x ■ x) = 1 is introduced to allow a Hamiltonian representation for the 
equations of motion (see Section 3.6 of m)- Then, the equations of motion are given 
by the Euler-Lagrange system with constraints, 


(3.4) 


^(avi(«.x,v)) - dx(t,^,v) - ASx9(x) = 0, 


where A is the Lagrange multiplier. After substitution, equation (3.2) simplifies to 

(3.5) crv(x ■ x) — ct(v ■ v)x — VyJJ{t, x) — 2Ax = 0, 

where A is unknown and we have used the fact that x ■ v = 0 for x G M^. To determine 
A, we take hrst the scalar product of x with the left hand side in (3.5) and obtain 

(3.6) cr(x ■ v)(x ■ x) — (v ■ v)(x ■ x) — X ■ x) — 2A(x ■ x) = 0. 

Since the particle is constrained to M^, we can differentiate equation (3.3) twice with 
respect to time to get 


V ■ V 


X-V = 0. 


Using this equation together with Euler’s formula (3.16), we can simplify (3.6) to 

-(v ■ v)(x ■ x) - (v ■ v)(x ■ x) - 2A(x ■ x) = 0, 

from which we get A = —(v ■ v). Substituting this value of A and x ■ x = a into (3.5) 
leads us to the equation 

V = VxU(x) — (j(v ■ v)x, 

a remark that completes the proof. □ 

The following proposition gives the equations of motion in local coordinates. 

Proposition 3.2. In local coordinates, the equations of motion for a particle with 
position q = + 6eg G and velocity p = oJaBa + snaujgeg G TM^, under the 

effect of a potential function U : M x —)■ M are given by the system 


(3.7) 


a = Ua, 

9 = UJg, 

^ uja = daU{t, a, 6) + cdg sn a csn a, 


uje = 


sn^a 


dgll(f, a, 6) — 2uaUJe ctn a. 


Proof. The Lagrangian for our system is given by the difference between kinetic and 
potential energies, i.e. 


(3.8) 


L{t, a, 6, a, 6) = -(d^ + ^^sn^a) + [/(t, a, 6). 


Substituting this into the Euler-Lagrange equations, 
(3.9) 

yields the desired equations. 


^(daL) - daL = 0 and ^(dgL) - dgL = 0, 


□ 
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The following obvious result expresses the equations of motion in the context of the 
Hamiltonian formalism. 

Remark 3.1. In extrinsic coordinates (x, v) G TM^, the equations of motion for a 
particle of mass 1 moving on can be written in Hamiltonian form as 

X = d^H 
V = -d^H, 


where if(t,x, v) = ■ v)(x ■ x) — f/(f, x) is the Hamiltonian function. 

3.2. Gravitational potential and gravitational field. To dehne the Vlasov-Poisson 
system on M^, we need to obtain a suitable solution representation to Poisson’s equation 
on and This representation is provided below in terms of the potential function 
U. 

Proposition 3.3. The gravitational potential function U at a point x G due to a 
spatial mass distribution p = p(t, x) is given by 


(3.10) 




M2 


p{t, y) log 


1 + (jx ■ y 
cr - X ■ y 


dy. 


Proof. The proof is trivial given that the fundamental solution of the Laplace-Beltrami 
operator on the manifold is given by (see for example 00) 


(3.11) 


C/(x) = — logctn 


<i(x,y) 


where d is dehned as in (2.2). In particular, it satishes 

6 {a - a') 0 6 {e - 6') 


(3.12) 


Ax5(x) = 


sn a 


where 5 is the usual Dirac delta distribution and ® denotes the tensor product. Thus 


(3.13) 


[/(x) = 


^(x)p(y)dy. 


Now, since (i(x, y) = csn ^(ax ■ y), and thanks to the trigonometric identity 


d 1 + csn d 

ctn- = a- --, 

2 1 — csn d 


the potential simplihes to 
(3.14) 
as desired. 


= 4- 


log 


m2 


1 + crx ■ y 
cr - X ■ y 


P{y)dy, 


□ 


Observe that if we extend the potential U as 

X \ _ 1 


(3.15) 


D(x) = U 


p(y) log 


/ ||x|| + crx ■ y 
Vcr||x|| - x-y 


dy, 
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then U becomes homogeneous of degree 0, and Euler’s formula yields the identity 
(3.16) X ■ x) = 0. 


The physical interpretation of (3.16) is that, in an extrinsic coordinate system having 
the origin at the centre of M^, the gravitational force acting on each particle is always 
orthogonal to its position vector x. Therefore we can conclude that if particles are 
initially located on with velocities tangent to the manifold, then they will remain 
on for all time. The proof of this fact, obtained in Section 3.7 of m for the discrete 
Newtonian case, can be easily extended to our problem. 

Proposition 3.4. The gravitational field present at position x G is given by 


(3.17) 


F{t,x) VxU(t,x) = 


27ra 


m2 


c^(x-y)x 

^p{y)dy- 


yy 


Proof. We begin with (3.10) and extend the expression by homogeneity so that it be¬ 
comes 


(3.18) 




M2 


p{t, y) log 


|x|| -|- crx ■ y 


O' X 


xy 


From here, we use the gradient as in (2.3) to find 

1 


(3.19) 


V,C/((,x) = 


2 'Ka 


|x||y 


O' X 


y)x 


m2 


y)^ 


dy. 


p{y)dy. 


Invoking ||x|| = 1 yields the required result. 


□ 


If we consider a specihc mass distribution, it is also possible to use Gauss’s law to 
calculate the gravitational held without hrst knowing the gravitational force function. 
In fact this method proves to be useful in the case of symmetric mass distributions, as 
shown in the following example. 

Example 3.1. Let a mass distribution on be given by 


p(a, 6 ) := p(t, a, 6 ) = 


h(a)®5(0) 


sm a 

i.e. a point mass located at the north pole of the sphere, as shown in Figure Gauss’s 
law in two dimensions says that the gravitational hux at any point g G is proportional 
to the mass enclosed by a Gaussian curve passing through the point, i.e. 


(3.20) 


F • n = m. 


'C 


is the normal vector to the curve G, 


where F denotes the gravitational held, n G TqS 
and m is the enclosed mass. 

If we choose a Gaussian curve C as pictured in Figure the spherical symmetry of 
the mass results in signihcant simplihcation of (3.20). Since the held generated by the 
mass intersects the Gaussian curve perpendicularly within and the magnitude of 
the held is the same at each point on the curve, our equation reduces to 

27rsinQ;|F| = 1. 
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Figure 1. Illustration of a point mass located at the north pole of 



Figure 2. A Gaussian curve, G, chosen for Example 3.1 


Solving for F yields 

(3.21) F = -e„. 

271 sm a 

Alternatively, we can calculate the same held by taking the gradient of the gravitational 


force function for the point mass. To do so, we hrst hnd the force function using (3.10) 
with ([^, 


r‘27r /‘TT 


U{t,a,e) = - 


'0 JO 


5{a') ® 5{e') 

sin a' 


log 


= 4ii'°S 




1 + x(q,6>) •y(0,0) 
1 - x(a,6') ■ y(0,0) 


1 + x(a, 9) ■ y(a', 9') 
_1 -x(a!,6') ■y(a',6'') 

1 , /1 + cos a 

= log ^- 

471 V 1 — cos a 


sin a'da'd9' 


I ct 
cotl- 
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then use ( 2 . 8 ) to take the gradient and obtain 


V^U{t,a,e) = 

Zn 


a , 

cot — e„ = 


27r sin a 




which agrees with (3.21). 


3.3. The Vlasov-Poisson system on spaces of constant curvature. According 
to kinetic theory and Liouville’s theorem (see, e.g., |22]), the equation that governs 
the motion of a continuous particle distribution with no collisions is given in local 
coordinates by 


^f{t,a,e,Uc,,ug) = 0 , 

at 


where / = f{t,a,9,Ua,uJe) ■= fit,^{a,9),\(ua,uje)) and / is the phase-space distribu¬ 
tion function in extrinsic coordinates. Using the chain rule, this equation becomes 


dtf + adaf + 9 def + + ug d^J = 0. 


Employing the equations of motion (3.7), we can write the above equation as 


(3.22) 


0 = dtf + uJadaf + (^edgf + (9„f/ + oi^sn acsn a)d^^f 

dgil - 2uau:ectna \ 


sn^a 


which we call the Vlasov equation on or the curved Vlasov equation. Using the 
same approach in extrinsic coordinates yields the equation 


(3.23) dtf{t, X, v) -h V ■ Vx/(t, X, v) -F [VxU(t, x) - cr(v ■ v)x] ■ Vv/(t, x, v) = 0, 

where Vv/ = d^^f ei + d^^f e 2 -l- crdv^f 63 , and we require that x ■ x = a and x ■ v = 0 , 
such that the particles remain on during the motion. When we couple the Vlasov 
equation to our solution representation of Poisson’s equation and add a compatibility 
condition between our spatial density and phase-space density, the result is a closed 
system, which we call the gravitational Vlasov-Poisson system in spaces of constant 
curvature, or for short the curved gravitational Vlasov-Poisson system. 


(3.24) 


r+ V ■ Vx/ + (Vxf/ - a(v ■ v)x) ■ Vv/ = 0, 




p(t,x) = 


/Tx(M2) 


p(Uy)iog 

M2 

/(t,x, v)dv. 


1 -|- crx ■ y 
cr - X ■ y 


dy, 
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In local coordinates, system (3.24) takes the form 


(3.25) 


'dtf + Uadaf + Ued0f + {dJJ + Wg snacsna)c}^^/ 

-i^deU - 2a;„a;ectn a ] d^J = 0 

J 

1 + crx(Q!, 6 ) ■ y («', 9') 

pyi, a , (7 j log 
m2 




^sn 

p{t, a', 9') log 


p{t,a,9) = 


/Tx(M2) 


a — x(a;, 9) ■ y{a', 9') 
f(t, a, 9) sn a duadue- 


sn a' dad9 


Remark 3.1. In [25] an approach in the spirit of geometric mechanics was used to 
define the Vlasov equation. In this sense, we can first define the Hamiltonian as 

H{q,p) = ]^\p\^+ U{q). 


Then, the conservation of a distribution function f(t,q,p) along the phase-space trajec¬ 
tories of H, i.e. the law ^f = 0, can be written, using the canonical Poisson bracket 
as 

dtf = {f,H}, 


which gives the Vlasov-Poisson system if we require that the potential U solves Poisson’s 
equation. 


-AM 2 U{q) = p{q) : = 



f{t,q,p) dp. 


Although this form of the Vlasov-Poisson system in curved spaces is elegant and quite 
natural, its analysis seems to be difficult. For this reason, we will focus in the sequel 
on a special conhguration of the mass distribution. 


4. Initial data along a geodesic 


In the rest of this paper we focus on the 1-dimensional case, and assume an initial 
distribution on in which the particles lie initially on a geodesic, which is a great 
circle in and a hyperbolic great circle in In our model of hyperbolic geometry, 
a hyperbolic great circle is a hyperbola obtained by intersecting the upper sheet of 
the hyperboloid of two sheets with a plane through the origin of the extrinsic coordi¬ 
nate system. Therefore, we study a conhguration of particles that obeys the following 
conditions: 


1. the particles move on the geodesic G, where 


G : = 


{{Xi,X 2 , Xs) \ xj + xl = 1, Xs = 0} , 
{{xi,X 2 , Xs) \xl-xl = -1, = 0} , 


for 

for m2 


§2 

e2. 


2. the velocity of each particle is always in the eg-direction for M2 = §2 
eQ,-direction for M2 = ]Hl2, 

where the geodesic G has been chosen for convenience and without loss of generality. 
Notice that, since the gravitational force on each particle is directed along the geodesic. 
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if the particles are initially aligned on G with initial velocities along that geodesic, then 
they remain on G for all time. 

For convenience, we consider 


r [0,27r), if M2 = §2 

I {—oo,+oo), ifM2=]Hl2, 

and for any functions (f){x,oj),'ip{x,u)) with x G 1 , 0 ; G M, we dehne 


<(t),'ip> = 



da duj. 


Let us further derive the equations when the motion of the particles is restricted to 
the geodesic G. 


4.1. Motion on a great circle of §2. Let us notice that in the spherical case we have 

G = {x(a, 0) : a = 7r/2} . 

This allows us to take the phase distribution after restriction in the form 

(4.1) f{t,a,e,Uo,,U0) = 0 S{ua)g{t,e,U0), 

sm [a) 

where is a distribution on G. If we let pg(t,9) = J^g(t,9,u;0)du0, then the following 
proposition, which is our hrst notable result, provides the equations of motion restricted 
to G. 


Proposition 4.1. If the phase distribution is taken as in (4.1), then g and pg satisfy 

dtg + ued0g + F{t, 9)d^gg = 0, 

F{t,9) = W*0d0Pg, 


(4.2) 


pg{t,9) = / g{t,9,U0)du0 


in the distributional sense, where W{9) = ^ log | cot(|)|, 9 G [0,27r), and *0 denotes 
convolution in the 9 variable. 


Proof. Let 0 G C“([0,27r) x M) be a test function. Then there exists <F G C“([0,7r] x 
[0, 27r) X M X M) such that = 4>- Noticing that, after restriction, oIq, = 0 and 


daU\a=^ = 0, substituting / into (3.25) and testing with $, we have for the hrst term 



2 ?! /* 7 r 


271 


/ dtf^dad9dujadu0 = dtg4>d9duj0. 

'0 Jo Jr Jo 

The computations of the other terms are similar and in the end we obtain 
(4.3) < dtg + oj0d0g + deUgd^^g, 0 >= 0, 

where Ug = U\a=-^. 

We need to obtain now the explicit form of Ug. Since for x(q;, 6*), y(a', 9') G § 2 ^ 
X ■ y = sin a sin a' cos(6* — 9') + cos a cos a', 
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we have that 


U{t,a,e) 

n , ni\'ci I / \ , 1 + sinctsinct'cosf^ — 0') + cosacosct' , , 

/ pg{t,0)o{a —71/2) log ^^-:- 77 ^- 77 ^-- :da dO 


1 

471 


0 ^0 

1 / n/M 1 + sinct cosf^ — 6^') 

= ^ / Pgit, e')\og-^ —^- \,-^de'. 


1 — sin a sin a' cos(6' — 6 ') + cos a cos a! 


4vr _ 
Therefore, 


0 


1 — sin a cos(6* — 6*') 


Ug = U{t,7l/2,d) 

1 r" / o/M 1 + cos(e - 6') 
— / Pg{t, 9') log-. - 


r-27r 


— / Pg{t, 9') log 


1 — cos(6' — 9') 
{9-9' 


If we dehne the kernel 


we can write that 


IT(0) = —log 


cot 


9 

cot - 
2 


d9'. 


□ 


d0Ug = W *e dePg. 

Setting F{t, 9) = delJg completes the proof. 

4.2. Restriction to a hyperbolic great circle of Let us note that, in the 
hyperbolic case, 

G = {x(a, 9) : 9 = 7r/2, a G {—oo, +C)o)} . 

We can therefore take the phase distribution after restriction in the form 


(4.4) 


f{t,a, 9 ,Ua,u} 0 ) = ® S{ue)g{t,a,Ua), 


sm^{a) 

where g is a. distribution on G. If we let Pg{t, a) = g{t, a, Ua)dUa, then the proposition 
below provides the equations restricted to G. 


Proposition 4.2. If the phase distribution is taken as in (4.4), then g and pg satisfy 

' dtg + ojadag + F{t, a)duj^g = 0, 

F{t, a) = W *a dapg, 

pg{t,a)= / g{t,a,Uc:,)dUa 

Jk 

in the distributional sense, where W{a) = ^ log | coth(|)|, a G M, and *a denotes 
convolution in the a variable. 

Proof. As in the spherical case, let 0 G x M) be a test function. There exists 

<h G (^“(M X [0, 27r) X M X M) such that = <p. Noticing that after restriction, 

U 0 = 0 and d 0 U\ 0 =^ = 0, substituting / into (3.25) and testing with $, we have 
(4.5) < dtg + ujadag + dJJgd^^g, 0 >= 0, 
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where Ug = 

Since for x(q!, 9),y(a', 6') G 

X ■ y = sinh a sinh a' cos{9 — 9') — cos a cos a', 

we can write that 
U{t,a,9) 

2 r /•27r 

4?! 



Pg{t,a)6{9' - 7r/2) log 


1 — sinh a sinh a' cos{9 — 9') + cosh a cosh a' 
— 1 — sinh a sinh a' cos{9 — 9') + cosh a cosh a' 




1 — sinh a sinh a' cos{9 — 7r/2) + cosh a cosh a' 
— 1 — sinh a sinh a' cos{9 — 7r/2) + cosh a cosh a’ 


da'. 


Therefore, 




Ug{t,a) = U{t,a, -) 


- //»(*.«') log 
i /p»(‘.o')log 


1 — sinh a sinh a' + cosh a cosh a' 
-1 — sinh a sinh a' + cosh a cosh a' 
1 + cosh(Q; — a') 


da' 


T + cosh(Q! — a') 


da' 


Pg(t, a') log ^coth^ 
coth 


1 

4?! _ _ 

27r ./» 


a — (X 


a — a 


da' 


da'. 


Setting the kernel 


and the force 


completes the proof. 


lT(a) = —log 


1 “ 

coth — 


F(f, a) = dJJg, 


□ 


We can combine the resnlts of Proposition 4.1| and Proposition |4.2| by replacing both 
{9,oj0) in Proposition 4.1 and {a,Ua) in Proposition 4.2 by {x,v) G / x M. Doing so 
prodnces the following proposition. 

Proposition 4.3. Let f be a phase space distribution along any geodesic G of W?. 
Then the system restricted to G becomes 

' dtf + vd^f + F(f, x)d^f = 0, 

Fit, x) = W * d^p, 


(4.6) 


p{t,x)= / fit,x,v)dv, 


where ix,v) with (x,n) G / x M, hP(a;) = ^ log |ctn (|)|, and * denotes the convolution 


m X. 
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Several quantities are conserved along the solutions of the above equations, and the 
next result summarizes them. 


Proposition 4.4. Let f{t,x,v) be a solution of (4.6) such that f has compact support 
in V, has compact support in x for and is 2n-periodic in x for Then 

system (4.6) can be alternatively written as 

(4.7) dtf+ vd^f + da^Ud^f = 


x-y 


dy and p{t,x) = / f{t,x,v)dv, and the 


where U{t,x) = — / p(t, y) log ctn 
271 

following quantities are conserved: 

(i) the total number of particles, 

N:= f{t,x,v)dxdv, 

Jr Ji 

(a) the total mechanical energy, 

f{t,x,v)v‘^dxdv— U{t, x)p{t,x) dx, 


E:=- 

2 



(Hi) the entropy. 


(iv) the Casimirs, 


S := 



C : = 


f{t,x,v)\og[f{t,x,v)] dxdv, 


A{f(t, X, v)) dxdv, 



where A is an arbitrary smooth function. 


Proof. To show that the total number N of particles is conserved, we integrate (4.7) 
over the phase space as follows. 


0 = 




[dtf + vd,,f + d:AJd^f] dx dv 


dtfdxdvP vdxfdxdvP d^Ud^fdxdv 




^ fdxdv-\- vdxfdxdv-\- / / d^Udyfdvdx 





jt I I 

— fdxdv, 



f dxV dx dv 


f dy^{dxU) dvdx 



where we used that / has compact support in v, that /(O) = /(27r) in the positive 
curvature case, and that / has compact support in x in the negative curvature case. 
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For the conservation of total mechanical energy, 

fv^dxdv — 


E = - 
2 



Updx, 


we multiply (4.7) by and integrate over phase space. 

Since / satishes relation (4.7), by the chain rule so does the Casimir A{f). Therefore, 
we can integrate 

dtA{f) + vd^A{f) + d^Ud^A{f) = 0 

over the phase space to yield our conservation law. Entropy is a Casimir, so its conser¬ 
vation follows directly from the conservation of Casimirs. □ 

We are now also in a position to prove the following result. 

Proposition 4.5. Any distribution of the form f{t,x,v) = f^{v) is a spatially homo¬ 
geneous equilibrium solution of equation (4.6). 


Proof. By dehnition, any stationary solution must satisfy (4.7) with dtf = 0, i.e. 

(4.8) vda:fF dyf = 0. 

Consider a spatially homogeneous distribution function / = f^{v). For this form of /, 
we get 

(4.9) = 0, 


so the hrst term in (4.8) is 0. We can use the expression of p in (4.6) to calculate 

(4.10) p{t,x) = J f^ dv = p^ (constant). 

Consequently, we obtain that the force due to the homogeneous distribution is 

(4.11) F = W*d^p^ = 0, 


since is a constant. Therefore the second term in (4.8) vanishes and we conclude 
that f{t,x,v) = /°(n) is a spatially homogeneous equilibrium (stationary) solution to 
equation (4.7). □ 

5 . Local well-posedness 


The purpose of this section is to construct a solution to (4.6) in the class of analytic 
functions such that to apply to it the Penrose condition for linear stability. In a forth¬ 
coming paper, we will study the existence of weak solutions and strong solutions local 
in time in the class of Sobolev spaces. This is in contrast with the Vlasov equation 
given by the Dirac potential studied in [19] , where the singularity is stronger. Since we 


are working with system (4.6), we rewrite it here in a slightly different but equivalent 
form. 


(5.1) 


^ dtf + vd^f -P{W * d^p)d^f = 0 
P= [ fdv. 


W{x) = — log 


ctn|| 
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We will construct our solution in the class of Gevrey functions G*. For s > 1, set 

G"(x) := {/(x) I Vk = (fci, k 2 ) e 3Mk > 0 s.t. \D^g\ < M^{k^\k 2 \y] . 

Note that is the set of analytic functions and C Let P be the cone of 
polynomial functions with positive coefficients. Fixing s > 1, we dehne the operator 

Da.,S p 

n\ 




{n — a)\ 


where a,n eN. 


Remark 5.1. When s = 1 the operator 0°“’^ reduces to the classic a-derivatives operator 
w.r.t \, which will help us define a suitable function space to find the solution. 

For any function f{t,x,v) and positive real numbers {Xq,K), dehne 

h., = d^Af, l/(i)lA:= 7T^I/«U-, l/(i)lA,a:=-D“'‘l/WU, 

k+l>0 ^ ■ '' 


X{t) = Xo-{l + K)t, Hf{t) -.= ^2 


\f\x{t),ait) £V / X __ 2 sr , i\s-l l/U(l))“(^) 


a>0 


(a!) 


2s 


W(t) 


a>0 


a! 


\2s 


Note that all the above norms still make sense when the functions depend only on the 
variable a: or n by setting A; = 0 or / = 0, respectively. 

Let us denote G^ = ||hF||Li and assume that the initial condition ffix^v) of (5.1) 
satishes 

(5.2) ho-.= Hffitf) < oo. 

For all M > ho > 0, we set M = In and for all T > 0, we dehne the space 

'^T,M = {f{t,x,v) : sup Hf{f) < M, [ Hf{t)dt < M}, 
ie[o,T] Jo 

which we endow with the norm 


T = sup Hf{t) + / Hf{t)dt. 
te[o,T] Jo 

We are now in a position to state our second notable result. 


Theorem 5.1. Let s > 1 and let fi G G^ be the initial condition of (5.1) satisfying 
ho < cxD. Then there exist a time T > 0 and a constant M > ho such that (5.1) has a 
solution f G 'Ht,m C G^. 

The proof of this theorem is tedious and goes through several steps. First, we need 
the following propositions and lemmas. 
















The Gravitaional Vlasov-Poisson System in Curved Spaces 


19 


5.1. Preliminary estimates. In this subsection we will set the background for the 


proof of Theorem 5.1 Our hrst result generalizes the Leibniz rule and provides some 
useful inequalities for the operator 

Proposition 5.1. The operator is linear and for any a,b eN, p(A),g(A) G P, 

(5.3) = D^+^’^p{X), 

a 

(5.4) (1 + < YXCtYD'‘-‘p(\)Dl-’‘q{X) < £>“’*(p(A),(A)). 

k=0 


Proof. The proof of (5.3) is straightforward. For (5.4), it is sufficient to prove that the 


inequalities hold for p{X) = A'", q{X) = A"'. For this, notice that 

a a 


k=0 


k=0 


\fc=o 

Meanwhile, since 

1 “ 


a—k \ ^m-\-n—a 


(m + n)! 
(m + 77 , — a)! 


X^+n-a ^ D-,^^p(^X)q{X)). 


o T 1 




k=0 


k=0 


we have 


^(CS*B‘'*p(A)B‘-‘'*,(A) = 5;(a!C‘Cr‘)‘A 


m+n—a 


k=0 


k=0 


> 


(a + 1)‘- O; a!C‘Cr‘ A”+"-“ = (a + l)‘-*C“'‘(p(A),(A)). 


, fc =0 


This remark completes the proof. □ 

It is important to note that estimates on the distribution function f(t,a,u) in Hf or 
Hf norms automatically yield the same estimates for the corresponding density function 
p. Indeed, since a{v) > 0 and J^a{v)dv < 1 , we have 

\d^P\L^ = I [ ad^fdv\L^ < I [ (^{v)dv\\d^f\L^^^ = • 

t/ M. M 

The following result enables us to estimate the distribution function when multiplied 
by the weight function 7 . 

Proposition 5.2. Assume 'y{v) satisfies 

(5.5) C^:=^Y^ C[X^^ sup (|7n|L-) < 00 , 

fc>0 0<l<k 0<n<k 

then for any f{t,x,v) we have that 

H.yf < C-yHf and H^f < C-yHj. 
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Proof. First, using the definition of the norms and the Leibniz rule, we have 


a>0 k-\-l>a 

= EE 


a>0 k-\-l>a 

^ EE 

a>0 k-\-l>a 


1 A^+'-“ 

( {k + l)l 

(a!)2* (A;!/!)* ' 

(A; — a)! 

1 

f {k + l)l 

(a!) 2 * {klliy ' 

1 

+ 

1 

f {k + l)l 

(a!) 2 * {klliy ' 

1 

+ 


|(7/)a:,;|l“„ 


6=0 


;,6 L ° 


0<6<Z 


sup (|7n|L“) 

0<n<l 

EEE 

0<b<l a>0 k+l>a 


1 Xk+b-a / ^k+b)\ 


{a\y^ {k\b\y \{k + b-a)\ ' 


C'f 


b\ (k + iy. {k + b-a)\ 
n{k + b)\ {k + l-a)\ 


y ^ sup (|7nU-) 


0<n<i 


< 


EEE 

0<b<l a>0 k-\-l>a 

sup (|7n|L-) 


1 Xk+b-a / (k+b)\ 


(a!)2^ {klbiy \{k + b — a)\ 


\fk,. 


6 L2 


0<n<l 


Second, since b < I, for each fixed I we have 

1 A^+^-“ / {k + b)\ 

a>0 k-\-l>a 


a\y^ {k\b\y \{k + b-a)\ l 


< 


EE 

a>0 k-\-b>a 


1 A^+'’-“ / {k + b)\ 


(a!)2* {klbiy \{k + b-a) 


\fi 


fc,6|L“ < Hf. 


Therefore 


H-,, < E E «/ ■ 7a‘ sup (|7„U«) < C,H,. 


l>0 0<b<l 

The proof for H.yf is similar. 


0<n<i 


□ 


The next result is the first step in the key estimate which shows that the contraction 
map is well-defined. 


Lemma 5.1. Given p{t,x) and x{v) satisfying (5.5), let g be the solution to 
(5.6) dtg + vd^g + {W * d^p){d^g + xg) = 0, 


with g{0,x,v) = fi{x,v). Then we have that 

dt\9\x,a < U^|5'U,a + Al^flA.a+l + (|p| a,i l^'|A ,1 + |p|A, 1 175'Ia)• 
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Proof. Differentiate (5.6) k times w.r.t. x and I times w.r.t v and notice that 


dt{gk,i) + vd^{gk,i) + (ID * d^p)dy{gk,i) 


k-l 


+ 5; ct(w »aj-”+V)(9„,,+i) + Y. 


m=0 


m=0 


Multiplying by ’ summing over k + l>a, and using the 

method of characteristics, we obtain the inequality 


(5.7) 


(fcK^-1)!)^ \{k + l-a)\ 


\gk+i,i- 


1 n; 


(5.8) +C^Y1 Y1 


k>a—l m<k—l 


{l\m\{k — m)\y 




(5.9) +C.,Y 


k>a—l m<k 


(l\m\{k — m)!) 


y 79 )™,iU.~jsT’”+vii, 


The hrst term on the right hand side can be controlled by changing fc + 1 to fc and I — 1 
to /, as follows, 


_ X^+l-a^^k + m^ 

Q {km{k +1 - 


V 2^((l/2)a + (l/2)(fc - a)y ^ \gki\L^ 


k>a—l 


< 2-v V ^'^'-^jjk + iy.y 

< ayg\x,a + Mg\x,a+i- 


+ 2- 5^ 


k>a—l+l 


x^+^-y{k + ly.y 

{km{k +1 - a - lyy 


■\g^ 


'k,l\L^ 


The second term is (5.8) = CwD^^’^Y, where 


i'' = E E 


k>a—l m<k 


^ (l\m\{k — m)!) 


l9«..i+ili?jaT”+Vli 
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Because y is a polynomial in A with positive coefficients, we can get an estimate on 
D‘^Y just by estimating Y itself, 


Y = 


< 


< 


E 

E 

k>a—l m<k—l 

E 

E 

m+/>a- 

-k k>2 

E 

E 


:l9™.i+ilirJSL"+VlL 


(l\m\{k — m)\y 

yk-\-l+m—2 


((/ - l)!m!(A: - 1)!)^ \ 9 mAL^Jd,p\Lo 


yl+m 1(^/ _|_ ffiY \ 


k-1 


m+l>a—k k>2 






-\dlp\ 


L2° 


yk-l . 

^ - {llmiy - 2^ JYP - - IpIa,i|£/|a, 


m+/>l 


k>l 


iik-iy.y 


Hence 

Similarly, 

where 


(5.8) = C^D'^'^Y < C^D‘^’y\p\x,i\9\x,i). 


Y = 


E E 


(1^ < C^D^’^Y, 

^ \{\9)mAL^M~^^" 


k>a—l m<k 


(l\m\{k — m)\ 


(^\l° 


\m+/ 

E E 


A 


k—m 


k>a—l m<k 

Consequently we have 


(mW.y 


I LZ 


\d^ < IpIa,i|7P|a- 


{{k — m)\y 


(5.9)<aD“’^(|pU,iN|A) 


The proof follows by combining estimates (5.7), (5.8), and (5.9). 


□ 


1 


In ^ are given. 


18 '>C™+ 2 C.y ho 


Lemma 5.2. Assume M > 0, 7 (u) satisfying (5.5) and M = 

Let \{t) = Ao — (1 + K)t, 0 < f < T. Then there exists K such that for any p with 


sup Hp{t) < M and 

0<t<T 


Hp{t)dt < M, 


the solution to (5.6) is such that g G 'Ht,m o-nd 

dtHg < (Ao - iC + (4 ■ 9* + CyC^Hp)Hg + (18^ + 2CyC^HgHp. 

Proof. By Lemma [13 

^t|5'|A(t),a < + •^(^)l5'|A(t),a+l + (|p| A(p,l [S'I A(t),l + |pU(p,l ItpI A(t)) 

— (1 + iL)|5f|A(i),a+l 

= ^y9\gx(t),a + (-^0 - 1 - K)\g\x{t),a +1 + C'^-D“’^(|p|A(p,i|5'|A(i),i + IpIa.iIapIa)- 
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Multiplying the above inequality by and summing w.r.t a yields 


( 6 . 10 ) 

(5.11) 

(5.12) 


< E q^lsUcu + (Ao - 1 - if) E 


a>l 


\r (a + 1 )* I I I 

+^w 2_^ 2 ^ -- \9\\{t),k+i\P\Ht),a-k+i 


a>0 0<k<a 


(a!) 




(a+irHc^a^r, , ,, 

-- \ig\\(t),k\P\Kt)A-k+i- 


a>0 0<fc<a 


(a!)^ 


It is easy to see that 


(5.10) <Hg + {\o-l- K)Hg < (Ao - K)Hg. 


Then changing the index k + 1 to k, we obtain 


= CvE E 


a>0l<fc<a+l 


(a!) 




a>0l<fc<a-2 


+Cto|p|A(t),l ^ 
a>0 


(a!)2- 

(a + 1)^ ^|fi'|A(t),a+l 


(a!) 


2s 


+ Cw\p\\{t),2 ^ 


(a + 1)"^ 


a>0 


(a!) 


2s 




a>l 


(a!) 


Finally, changing the index a — k + 2 to a, we obtain 


(5.11) = 


(5.13) 


CvEE 

a>4 k>l 


{a + k — 1) 


S—1 


{{a + k - 2)\y {{k - i)\{a - ly.y 


■ |fi'|A,A:|p|A,a 


„ (® T 1)* ^11 I I I I 

+ C^Pa.I 2^ t ,X2S I^U,a+l + h^|P|A,2 2^ |^|a, 

a>0 ^ ^ P 


a>0 

re 1 /it I I I I 

^ ^ ((«-l)0"" ^^^^’“■^2*(a + l)(a-l)^a2« 
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and 


(| 5 l^ = EE 


|^|A,fc 2 s IdU.a [ kk\{a-l)\ 


k>l a>4 


{a + k — 1) 


s-l 


{k\f^ (a!)2^ V(« + ^-2)!, 

k 1 CL 4 

^ \9\\,k/^ , ^ ^s_i^ 2 s lp|A,a ( k{k + l)\{a - l)\ 


(k\y 

k>l a>4 ^ ^ 

For fc > 3, a > 4, we have that 

A:(A: + l)!(a-l)! (a - 1)! 


2s 


V {a + k-2)\ 


■ 1 ■ 2 ■ 3 ■ 


4---{k + l)k 


< 6 . 


(a + A: — 2)! (a — 1)! a ■ {a + 1) ■ ■ ■ {a + k — 2) 

Similarly, for /c = 2, a > 4, we obtain that < 18, while for A; = 1, a > 4, 

are led to the inequality < 2. Then 


we 


(|5.13|) < 18'^C'^ifgFAp. 


Observing that for all s > 1 and a > 2 we have 

a(a + l)" ^ (a + l)*+i 


(5.15) 


< — 


1 / fl T 1 


2*(a + l)(a - l)^a2'< “ 2*(a + l)(a - l)^a2'< “ 2* V 
we can conclude that 


< 1 , 


(5.11) < 18^C^HgHp + AC^9^HpHg. 


Notice further that 




a>0 l<A:<a—1 


(a!)2« 

\s-l 


^ I I II ^11 I I 

+Cw\'jg\\ 2_^ / j\2s lpU,a+i + Cw\p\\A ^ (EnE IfpIa.o 

{^a^\Y-^{ClY. , , , 

-(y;Ti2i-l75'|A(t),fc|P|A(t),a-fc+l 

a>0l<fc<a-l 




Changing a — A; + 1 to A:, we are led to 


(5.16) 


(|^ < EE 


(a + k) 


S—1 


{{a + k — l)!A:!(a — 1)!)' 


■IbPlA.fclplA,^ 


a>2 k>l 

+CUHxgHp + HxgHg). 
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Meanwhile, 


|p|A,a [ A;!(a-1)! 


a>2 k>l 


((a - 1)!)2* (fc!)2- 


(a + A; — 1)! 




a>2 k>l 


((a-1)!)2* (A;!)2^ \{a + k - l)\ 


< 


a>2 k>l 

For a >2 and A; > 1, we have 

kk\{a — 1)! (a — 1)! 


^ I -|Ns-i |p|A,a \ig\\,k f kk\{a-l)\ \' 

{k'r [{a + k-iyj 


■ 1 ■ 


2 - ■ ■ k ■ k 


thus 


(a + A; — 1)! (a — 1)! a ■ (a + 1) ■ ■ ■ (a + A; — 1) 

_1 IpI A,a ITIAI A,fc 

((a - 1)!)2- (A;!) 


< 1 , 


( ra < 

a>2 k>l 
< CniHxgHp. 

We therefore obtain the estimate 


2s 


( |5-12[ ) < 2CwHxgHp + Cu,HxgHp < 2CwC^HgHp + CwC^HgHp. 
Estimates (5.10) and (5.11) together with ( 5.12[ ) imply that 

dtHg < (Ao - iF + (4 ■ 9* + C^)C^Hp)Hg + (18^ + 2C^)C^HgHp. 
For iF > Ao + (4 ■ 9® + C^)CwM , we have the estimate 

Hg{t) < = M. 

We can thus conclude that 

(/F - Ao - (4 ■ 9" + C^)M)Hg < (18" + 2C^)C^MHp - dtHg. 
Integrating from 0 to T yields 


cT 


(/C - Ao - (4 ■ 9" + C'^)C'^M) / Hgdt < {18^ + 2C.y)C^MM - Hg{T) + Hg{0) 


Thus for 


K > 


< (18" + 2C'^)C'^MM + 2M. 
(18" + 2C^)Cu,MM + 2M. 


M 


T -^0 T (4 ■ 9" + C'-y)C*^uiff, 


we obtain that 


Hgdt < M. 


We have thus shown that g G 'Ht,m and that the desired T, K satisfy 

^ ^ (18- + 2C,)C^MM + 2M. C.,)C.,M, 0<T< 


M 


1 +iF 
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This remark completes the proof. 


□ 


5.1 


5.2. Proof of Theorem 

that 7 (v) = a'{v)/a{v) satisfies 


Recall that a{y) >0 is chosen such that adv < 1 and 


fc>0 0<l<k 0<n<k 

Remark 5.2. Notice that functions a like above do exist. Indeed, take for instance 
a{v) = Ce~'^^, where 0 <C < \pk. 

Given a distribution function f(t,x,v), set p(t,x) = /go/dn and dehne a map by 
^{f)(t,x,v) := g{t,x,v), where g is the solution to 


dtg + vd,,g + {W * d^p){d.,g + 751 ) = 0 


with g{f),x,v) = fi{x,v). We will hrst show that there exist M,T such that $ maps 
T-Ltm into itself and then prove that $ is a contraction on 'Ht,m to finish the proof. 

First, observe that, thanks to Lemma 5T and Lemma |5.2[ the map $ is well de- 
hned. Next, let f{t,x,v),f{t,x,v) G 'Ht,m with f{0,x,v) = f{0,x,v) = fi{x,v), 
p = f^fdv,p = Jjg fdv and g = <F(/),^ = $(/). It is easy to see g — g satishes 

dt{g - g) + vd,cig - g) + (w * d,,p){d^{g - g) + jig - g) + {W * d,cip - p)){dvg + ig) = 0 . 


By Lemma [ 5 .1[ 


■^\g-g\\a < a"\g - g\x,a + {Xo - I - K)\g - g\x,a+i 

+Cu,D°''^{\p\x,i\g — g\x,i + IpUphig — g)\'r) 
+Cu,D°'’^{\p — p|a,i|^|a,i + \p — p|A,i|7dl7)- 


Estimates of (5.10), (5.11), and (5.12) of Lemma 5.2 imply that 


dtHg.g < (Ao - K)Hg_g R (18^ R 2C^)C^H,_gHp + (4 ■ r + C^)HpHg_g 
+ (18^ + 2C,)C^HgHp_p + (4 ■ r + C,)Hp_pHg 
= (18^ + 2C.y)C^Hg.gHp + (Ao - iF + (4 ■ O'* + C.y)Hp)Hg.g 
+ (18^ + 2C-,)C^HgHp.p + (4 ■ 9^ + C^)Hp.pHg 
< (18^ + 2C.y)C,,Hg_gHf + {X, - R N {4 ■ 9^ N C,)Hf)Hg_g 
+ (18* + 2C^)C,,HgH^_j + (4 ■ 9* + C^)Hj_jHg. 


If K 


/F — Aq — (4 ■ 9* + C^)M, we can write that 


|'g-/JC,,( 18 -+ 2 C+++s)d. 


dt 


< G^M(18* + 2C^)Hf_j + (4 • 9* + C.,)H^_^Hg. 











The Gravitaional Vlasov-Poisson System in Curved Spaces 


27 


Integrating over time, for 0 < t < T, yields 


^(A'/ 


< C„M(18‘ + 2C,) / Hj_i[s)d3 + (4 • 9* + C,)M sup 


0<s<t 


Consequently 

^ min (1, - ^11* < max (C^M(18^ + 2C^), (4 ■ 9^ + C,)M)\\f - /||,. 

Thus, taking t = T, we obtain 




M max((C'^M(18^ + 2C^), (4 ■ 9* + C^)M)) ~ 

ho min((l,ir-Ao-(4-9* +C'^)M)) 


Lemma 5.2 


Finally, in order to make <F a contraction, K,M,T should satisfy the conditions in 

Ao 


(5.17) K > ^ ‘2C^)Cy,MM + 2M ^ ^ 0<T< 

M 7/ > - - 


1 + 77 ’ 


(5.18) 


M max((C'^M(18" + 2C^), (4 ■ 9* + C^)M)) 
Yo min((l, 77 - Ao - (4 ■ 9* + C^)M)) 


These remarks complete the proof. 


□ 


6. Linear Stability 


In this section we will derive some Penrose-type conditions for linear stability around 
homogeneous solutions in the sense of Landau damping. Let us start by recalling from 
Proposition 4^ that the Vlasov-Poisson system restricted to a geodesic of is given 
by the equations 


' dtf + vd^f + F{t, x)dyf = 0, 

^ F{t,x) = W*d^p, 

p{t,x)= / f{t,x,v)dv, 

t. i/ M 

with {x,v) G 7m2 X M and potential W{x) = ^ log |ctn (|)|, where * denotes the 
convolution in x. 

We want to find conditions for linear stability around the homogeneous solution /°(n). 
Denoting by h(f, x, v) the fluctuation of / about /°, we can assume in the linearization 
process that the nonlinear term is negligible. In doing so, the fluctuation h becomes 
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the solution to the system 


( 6 . 2 ) 


( dh dh ,df 

F = W * 

ox 

p{t,x) = / h{t,9,v)dv. 


In the periodic case, the linear stability around the homogeneous solution /°(u) means 
that both the density and the force corresponding to the solutions of (6.2) converge 
exponentially fast to the space-mean of the density or to zero, respectively. More 
precisely, we have the following result in the case of 


Theorem 6.1. Assume that, in 

and the initial perturbation 

ho = ho{x,v) are analytic functions; 

- (/°)'(u) = 0(l/|u|) for large enough values of |u|; 

- the Penrose stability condition takes place, i. e. 


the stationary solution = f{v) of system (6.2 


(6.3) 


if 00 is such that (fA^oo) = 0, then p.v. 


iniv) 

V — 00 


dv > —1. 


Then there exist positive constants 6 and C, depending on the initial data, such that for 
t > 1 we have 


p{t,x) - 



2tt 


ho(x, v)dxdv 


<Ce and ||F(f,a;)||c’'(o, 27 r) < C'e 


-St 


C’’( 0 , 27 r) 


where ||M||c'-(o, 27 r) := maxo<n<r,o<x< 27 r \d^u{x)\ and r e N+. 


To state a similar linear stability result in the case of a hyperbolic circle in we 
need to redehne the norm of analytic functions in terms of the Fourier transform, see 
for example [26]. For any function /, dehne this norm by 


II/IUa := [ e«l|/(0|rf{, 

Jr 


where / stands for the Fourier transform of /. With this preparation, we can now state 
our linear stability result in 


Theorem 6.2. Assume that, in 

and the initial perturbation 

ho = hQ{x,v) are analytic functions; 

- (/°)'(u) = 0(l/|u|) for large enough values of |u|; 

- the Penrose stability condition takes place, i.e. 


- the stationary solution /° = /°(u) of system (6.2 


(6.4) for every cd G M, (f^Ylu) = 0 implies that 


p.v. 


'OO 
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Then there exist constants A', C > 0, which depend on the initial data, such that for t 
large we have 

2C C 

and \\F{t,-)y,,<-. 


The proofs of Theorem 6T and Theorem |6.2| rely on the following lemma about the 
decay of solutions of Volterra-type equations (see [3T] for its proof). 

Lemma 6.1. Assume that 0 solves the equation 

(j)(t) = a(t) + f K{t — T)4>{T)dT, 

Jo 

and there are constants c, a, Co, A, Aq, A > 0, such that the function a and the kernel K 
satisfy the conditions 

(i) \K{t)\<Coe-^o^; 

(ii) |A:^(0 -1\>C forO < Re^ < A; 

(hi) |a(f)| < ae~^^, 

where K^{^) is the complex Laplace transform defined by 

poo 

= / e^^K{t)dt for ^ E C. 

Jo 

Then for any positive A' < min(A, Aq, A), we have the inequality 

W«)|<Ce-y 

with 

C = O' A 


Coa 


2\/(Ao — A')(A — A') 

6.1. Proof of Theorem 6.1| , First, observe that the conservation of mass for the 
solution h of (6.2) is equivalent to the conservation of the zero mode of the density 
function p, 

p(f,0) = ho(0,0). 

In the following, we estimate each mode p{t, k) hrst, and then use it to show the 
convergence of pit, x) and F{t, x) with the exponential rate. The Penrose type stability 


condition comes naturally to satisfy the second hypothesis of Lemma 6.1 


In order to deal with higher modes, we hrst solve system (6.2) using the method of 
characteristics and a Duhamel-type formula. Setting S(t,x,v) := F(t,x)dvf^{v), we 
have 

(6.5) h{t, x,v) = ho{x — vt,v) — / S{t,x — v{t — T),v)dT. 

Jo 

Taking the Fourier transform h of h both in x and v gives 


( 6 . 6 ) 


h{t,k,rf) = ho{k,T] + kt)— / S{t, k,T] + k{t — T))dT. 

Jo 
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With the help of the identity 


we calculate that 


+00 

o = E(- 

.l)n+l 

cos( 2 na;) 

n 

n=l 


W{k) = \ 

[ M’ 

k odd 

1 

lo> 

k even. 


log 2 , 


Now, since S has the structure of separated variables, we have 

0 , 


S(T, k, „) = f(t, k)a,nv) = k)f („) = I ^ 


k even 


Substituting this expression into ( 6 . 6 ), we obtain 
(6.7) 

ho{k, 7] + kt), 


k even 


h{t,k,r]) = 


k 


ho{k,r] + kt)+ I + k{t - T))p{T,k)f°{ri + k{t - T))dT, k odd. 


Finally, the choice t] = 0 gives 
(6.8) 
where 


p{t,k) = ho{k,kt)+ / K{t — t)p{t, k)dT, 

Jo 

\ \k\tf^{kt), k odd. 


Clearly, for any given fc 7 ^ 0, ( 6 . 8 ) is a Volterra type equation. Since /° and ho are 
analytic, for large t we have 

\K{t,k)\ = and \ho{k,kt)\ = 


Assuming that the last hypothesis is also satished. Lemma 6.1 implies that 

with A' < min{Ao,Ai}. Now take r G M’*' and f > 1. Then there exists a positive 
constant C, which depends only on the initial data, such that 


P{t, k) - 



2-k 


ho{x, v)dxdv I 


k^O 


Choosing a constant 0 < h < A' such that for all h 7 ^ 0 and r G N"*" we have 

g-(A'-5)|fc|t < |h|-^-2, 

we can obtain the estimate 

CC 


p(C k) - 



2-k 


ho{x,v)dxdv I 


fc/O 


-5t 


2 











The Gravitaional Vlasov-Poisson System in Curved Spaces 


31 


where C* = ^2^=1 ^ Similarly for the force, we have 

d'.FM = Y. (■*)’■ vd". 


k odd 


and therefore, for f > 1, we have 




cc* 


^-st 


as desired. Finally, it only remains to check the second condition of Lemma 6.1 It 
basically means that when Re^ is located in a positive neighbonrhood of 0, the Laplace 
transform of the kernel K^{^, k) should stay away from 1. 

When k is even, K^{^,k) = 0, so it’s always away from 1. When k is odd, we 
evaluated the Laplace transform of K in time at .^ = (A — iu)k and obtained 


K^{^,k) = 








—iujkt tO 


r{v)dv 


(6.9) 


1 

W\ 


ifnv) 


-dv. 


iX + {v — oj) 

Moreover, if the stationary solution /° has the property that (f^Yiv) decays at least 
like (9(l/|n|), then (6 )^ implies that 


( 6 . 10 ) 




dv 


/]R v{iX - u + v) 

Let L{e) be upper half circle centered at 0 with radius £. Some simple computations 
show that 


1 

k 


V iX — u) + V 


dv = 


— lim 

k e-s-O 


^H-oo 


+ 


+ 


TTl 


k{iX — uj)~^ k £^o, 
(with V = 


JL{e) 

.0 ^ 


V iX — u + V 
1 


dv 


iX — oj + ee 

2711 




k{iX — u) ’ 


which means that (6.9) decays at least like 0(l/|a;|) as n —)■ cxd, uniformly for A G [0 , Aq], 
so it is enough to consider only the case when |n| is bounded. Hence, assume that 
|a;| < n. If (6.9) does not go to 1 when A —)■ O’*", by continuity there exists A > 0 such 
that (6.9) is away from 1 in the domain {|a;| < H, 0 < A < A}. Thus, we could only 
focus on the limit A —)■ 0"^. In order to compute this limit, we introduce the Plemelj 
formula (see IZSI). 


( 6 . 11 ) 


f(x) 


lim / 

y^o+ J_^ x-xo + iy 


—dx = p.v. 


r fix) 

’ —oo ^ ^0 


dx — i7if{xo). 
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Applying (6.11) to (6.9), we have that 
(6.12) lim K^iiX — iuj)k, k) = —-r^ ( p.v. 

A-5-0+ |fi:| 


foo j fO\i 


ifnv] 


V — LO 


dv - iTr{f ) (cd) 


We further need to hnd conditions such that (6.12) does not approach 1. First, if 
the imaginary part of (6.12) stays away from 0, then everything is hne. However, the 
imaginary part goes to 0 only if fc —)■ cx), in which case the real part also goes to 0, or if 
(/°)'(u) —)■ 0. So we only need to consider the case when (/°)'(a;) approaches 0. Hence, 
we need to require that 


(6.13) for every cd G M, (/^)'(cd) = 0 implies that — ( p.v. 

\k\ 


coo ( fO\i 


ifnv) 


V — u 


dv] ^ 1, 


a hypothesis that leads to the Penrose stability condition (6.3). 


□ 


6.2. Proof of Theorem 6.2[ Arguing as in the proof of the previous theorem, we just 
observe that hF(.^) = Atanh(^) (see [27]) and that 


(6.14) 


p{t,Cl = hQ{i,^t)+ / K{k-T)p{T,^)dT, 


where ^ tanh(^)/°(,^t). Moreover, the force is 




Computing the Laplace transform at C = (A — ico)^, we get 


K\C,0 = 


3 (A+iaj)^t j._1. / 


— tanh(^) / e y{v)dvdt 

^ ^ ./ro 


e(A+*‘^)«*-tanh(—) — 
2 ^ 2 


--tanh(y)(/ ) (u) 


(6.15) 




ifriv) 


iX + {v — cd) 


/ (/°)'(u)e-*«^‘dudt 

/r 

dv. 


Going though the same argument as we did in the previous section, we only need to 
consider the case when (/°)'(cd) approaches 0. Hence, we have to assume that 
(6.16) 

1 Pit / /■“ (f^y(v) \ 

for every cd G M, (/°)'((d) = 0 implies that — — tanh(y) ip.v. J - -dv ) ^ 1, 


V — u 


which is equivalent to 

(6.17) 

for every cd G M, (/°)'(cd) = 0 implies that ( p.v. 


.oo 


U — Cd 


dv ) 7 ^ 


2e 


tanh(^) 
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Since 




tanh(^) 


< —, we obtain the Penrose stability condition (6.4). Then the appli¬ 


cation of Lemma |6.1| provides us the exponential decay of the density 


for some constants C and A'. From this, we can conclude that 

llp(^c)ll^v = j <C f 

Jr Jr PJ 

and 

IT(*.')lhv= /h"'l«'tanh(^)|p(«,OI<i?<y / = §- 

Jr ^ ^ ^ Jr r 

a remark that completes the proof. □ 
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